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Abstract 

Magnetic behavior of a spin-1 Heisenberg dimer is analysed in dependence on both 
uniaxial single-ion anisotropy and XXZ exchange anisotropy in a zero- as well as 
non-zero longitudinal magnetic field. A complete set of eigenfunctions and eigen- 
values of the total Hamiltonian is presented together with an exact analytical ex- 
pression for the Gibbs free energy, longitudinal magnetization, longitudinal and 
transverse susceptibility. The obtained theoretical results are compared with the 
relevant experimental data of [Ni 2 (Medpt)2(^-ox)(H20)2](C104)2.2H 2 (Medpt = 
methyl-bis (3-aminopropyl) amine) . 
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1 Introduction 



One of the most challenging tasks emerging in the field of molecular mag- 
netism is to design molecule-based magnetic materials possessing predictable 
topology and desirable magnetic properties [1]. A successfully realized ratio- 
nal synthesis of low-dimensional molecular magnets would provide, among 
other matters, suitable model compounds for investigating magneto-structural 
correlations and in-depth understanding of cooperative quantum phenomena. 
Altogether, the simplest magnetic architecture among the molecular crystals 
have small magnetic clusters denoting the weakly-interacting assemblies of 
molecules formed by a few exchange-coupled paramagnetic centers [2]. 

It is of principal importance that a limited number of interacting sites in 
the discrete molecules often allows to model their physical properties within 
an exact quantum-mechanical approach that avoids a further approximation 
required for treating the extended solids. Among a large family of transition- 
metal coordination compounds are then dinuclear complexes the simplest model 
compounds which enable a precise description of quantum cooperativity and 
at the same time, provide an excellent testing ground for theories. Studies on 
their magnetic properties, moreover, played a crucial role in advancing the 
understanding of spin-exchange mechanism; it turned out that the nature and 
magnitude of exchange coupling depends on the ground-state electronic con- 
figuration of the constituting metal atoms as well as on the global topology 
and type of intervening non-magnetic atoms [3]. Besides, much effort directed 
towards the synthesis, design and analysis of dinuclear entities has been stim- 
ulated in connection with their wide application in various research fields. The 
metalloproteins with bimetallic cores often appear in the bioinorganic chem- 
istry [4], metalloenzymes (metalloproteins with a catalytic function) catalyse 
important chemical reactions at their metal centers in a variety of plants, 
fungi and bacterial species [5] and there is also a certain group of dinuclear 
complexes that shows an interesting photochemical behavior. 

In the present work, we will focus our attention on modeling the magnetic 
behavior of spin-1 dimer with a particular emphasis aimed at better under- 
standing of the influence of magnetic anisotropy. The main motivation to 
examine this model system can be related to a rich variety of existing homod- 
inuclear nickel(II) complexes, which can serve as suitable model compounds 
for the system under investigation. The most successful strategy used for the 
preparation of dinuclear complexes consists in an appropriate combination of 
the bridging and ancillary blocking ligands. The former ones are capable of si- 
multaneously binding two metal atoms in a close proximity and they also act 
as mediators of the superexchange interaction between the metal centers. By 
contrast, the latter ligands coordinatively saturate the metal centers in order 
to prevent the undesired self-assembly reaction providing the polymeric com- 
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plexes. Under these circumstances, the dinickel cores in neighboring molecules 
are well-separated from each other and there is only weak inter-molecular 
interaction through the hydrogen bonds and/or van der Waals contacts. 

From the magnetic point of view, the dinuclear complexes can be classified 
according to the nature of exchange coupling and the type of bridging ligand. 
At present, there exists a large variety of dimeric compounds in which the 
nickel centers interact through the antiferromagnetic (AF), or ferromagnetic 
(FO) superexchange coupling. In spite of a remarkable diversity of bridging 
ligands emerging in the dinuclear nickel-based complexes, the series of oxalato- 
[6] and 1,3- azido- bridged [7] compounds afford the generality of AF-coupled 
nickel dimers. On the other hand, almost all reported FO-coupled nickel dimers 
belong to the halide- or pseudohalide-bridged complexes [8] and among these, 
the coordination compounds with the 1,1- azido linkage provide their most 
extended sub-group [9]. It is worthy to notice that there also exists one rare 
exception of the FO dimeric compound, in which the dinickel core is embedded 
in the cavity of appropriate cryptand and the nickel atoms are linked via single 
azido group coordinated in the 1,3-bridging mode [10]. 

The organization of this paper is as follows. In Section 2, we will introduce 
the model system and we also recall the foundations for the occurrence of the 
magnetic anisotropy. This will be followed by a brief description of calcula- 
tion procedure: basic steps of derivation of the exact analytical relations for 
the Gibbs free energy, longitudinal magnetization, longitudinal and transverse 
susceptibility are listed along with a complete set of eigenfunctions and eigen- 
values of the total Hamiltonian. The most interesting analytical and numerical 
results for the zero- and finite-temperature quantities are presented and dis- 
cussed in Section 3. A subsequent part, Section 4, serves for a comparison 
of obtained theoretical results with the relevant experimental measurements 
performed on one representative example of the spin-1 dimeric compound. 
Finally, some concluding remarks are drawn in Section 5. 



2 Model system and its solution 

Before proceeding to an exact solution of the anisotropic spin-1 dimer, a few 
words should be addressed to a microscopic origin of the magnetic anisotropy. 
The global magnetic anisotropy of paramagnetic clusters is particularly as- 
sociated with the single-ion anisotropy and/or the g-factor anisotropy which 
are combined with more or less isotropic exchange coupling between the in- 
teracting metal centers [11]. It should be mentioned, however, that the main 
contribution to the magnetic anisotropy usually comes from the crystal field of 
ligands and hence, strongly anisotropic magnetic properties uprise just as the 
metal atoms reside in a crystal the sites of low (i.e. lower than cubic) symme- 
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try. Consequently, the spatial anisotropy in response to an external magnetic 
field is vigorously reinforced by: a) ligand asymmetry, when different types of 
ligands coordinate to each metal atom; b) geometric asymmetry, when there 
is an inequivalent geometric arrangement of ligands around each metal atom; 
c) coordination number asymmetry, when the different number of ligands are 
coordinated to each metal atom [12] . 

When the orbitally degenerate metal centers with a pronounced spin-orbit 
coupling build the metallic core of coordination compound, the exchange in- 
teraction between the magnetic ions becomes strongly anisotropic. This occurs 
for most of the rare-earth ions and also for some orbitally degenerate transi- 
tion metal ions such as octahedral Co 2+ and low-spin Fe 3+ cations. It is worth 
noticing that the exchange anisotropy differs from the single-ion anisotropy as 
well as the anisotropy of (/-factor. For example, the Yb 3+ - Cr 3+ exchange in- 
teraction in the YbCrBr 3 , - face-sharing bioctahedral dimer is well-described by 
the extremely anisotropic spin Hamiltonian: Ji = J xy (SfS% + SfS^) + JzS'^S^ 
with J xy =-4. 19cm" 1 and J2=5.16cm _1 despite both octahedrally coordinated 
Yb 3+ and Cr 3+ ions are by themselves magnetically isotropic [13]. Several re- 
cent studies have shown that an obvious exchange anisotropy can be found in 
the compounds with other orbitally degenerate metal ions as well [14]. 

On the other hand, the individual magnetic centers with an orbitally non- 
degenerate ground state usually necessitate the isotropic exchange interaction 
between them, since the anisotropic term becomes due to orbital quenching 
just of minor importance. In such a case, the origin of exchange anisotropy lies 
merely in the admixture of excited states through the spin-orbit coupling and 
thus, it can be viewed only as a small higher-order perturbation to the isotropic 
exchange [15]. The ground state of Ni 2+ ion in an octahedral environment 
is, for instance, orbitally non-degenerate and as such, it implies more or less 
isotropic intra-dimer interaction in the homodinuclear nickel-based compounds. 
Therefore, in an attempt to model the magnetic behavior of these compounds 
research has firstly focused on the theoretical interpretation based upon the 
assumption that the Ni 2+ ions are coupled by a purely isotropic interaction 
[16]. However, as it has been pointed out later on [17], the magnetic properties 
of octahedrally coordinated nickel compounds can be significantly affected 
by a rather large axial single-ion anisotropy. Bearing this in mind, Ginsberg 
et al. [18] took into consideration the uniaxial zero-field splitting parameter 
regarding the influence of single-ion anisotropy as well. 

Although the exchange anisotropy is in generality negligible for the most of 
homodinuclear nickel(II) complexes, we will incorporate it here into the effec- 
tive spin Hamiltonian to enable a deeper insight into the distinction between 
the effect of single-ion anisotropy and the effect of exchange anisotropy. To 
the best of our knowledge, such a comparative study has not been reported 
in the literature hitherto. Now, let us write the total effective Hamiltonian of 
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spin-1 dimer in the form: H = / H exc h + / H az f s + / H zee t, which is composed of the 
exchange term, axial zero-field splitting term and Zeeman term, respectively, 



7~(-exch — Jxy(SiS% + SfS^) + JzS^S^, (1) 

na,f. = D[(stf + (stf], (2) 

n zeet = -g z ^B(St + S z 2 ). (3) 

Above, S" and S% (a = x, y, or z) denote the spatial components of the local 
spin-1 operator on the metal centers 1 and 2, J xy and J z are the exchange 
couplings along the x, y and z directions, respectively, D labels the axial zero- 
field splitting parameter and other quantities have an usual meaning: B is the 
external magnetic field applied along the z-axis, /xb stands for Bohr magneton 
and g a (a = x, y, z) denotes the spatial component of g-factor. Notice that 
the negative (positive) sign of the zero-field splitting parameter D corresponds 
to an easy-axis (easy-plane) single-ion anisotropy. For easy reference, we will 
further refer to the z axis as to the principal axis. 

After performing elementary calculations in the usual matrix representation 
with a standard basis of functions (Sf = ±1,0 and S% = ±1,0), it is 

quite straightforward to find a complete set of eigenvalues and eigenf unctions 
of the total Hamiltonian: 



Ao,o = ~\j z + D-R, |* 0l o) = \[A+(\1, -1) + |1, -1)) - V2A_\0, 0)]; 
Ai )0 = —J z + 2D, |* 1>0 ) = -^(|1, -1) - | - 1, 1)); 

Ai,±i = -J xy + D T H, |* 1)±1 ) = -L(| ± 1, 0) - |0, ±1)); 



A 2 ,o = -\j z + D + R, |^ 2i(i 



i[A_(|l,-l) + |l,-l)) + v / 2A + |0,0)]; 



2 ' 1 ' ' 2 

A 2 ,±i = J xy + D T H, |* 2)±1 ) = -L(| ± 1, 0) ± |0, ±1)); 

A 2 ,± 2 = J Z + 2D T 2H, |* 2i±2 ) = I ± 1, ±1). (4) 

Here, we have introduced the expressions H, R and A± in order to write the 
relevant eigenvalues and eigenvectors in a more abbreviated and elegant form: 



H = g zl * B B; R=^(±-D) 2 + 2Jl y A ± = ^ R±{ \ ° l . (5) 

From here onward, one can readily proceed to the evaluation of basic ther- 
modynamic quantities within the standard thermodynamical-statistical ap- 
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proach. Substituting the eigenvalues (4) to a statistical definition of the parti- 
tion function leads, for instance, to the following exact analytical expression: 



Z = ^exp(-/3A Q ) = exp[(3{J z - 2D)} + 2exp[-[3{J z + 2D)} cosh(2/3#) 

a 

+4exp{-(3D) cosh{(3H) cosh(f3J xy ) + 2exp[(3{J z - 2D)/2] cosh(f3R), (6) 

with (3 = l/k B T. It should be stressed that the previous equation provides, 
because of the relation Q = — In Z, an exact solution also for the Gibbs free 
energy Q. Hence, some other important quantities such as magnetization, pair 
correlation functions and quadrupolar moment can be in turn calculated by 
differentiating the Gibbs free energy with respect to the relevant interaction 
parameter included in the total Hamiltonian. Thus, one easily obtains the 
magnetization per one site m z = —\ d ®jf H y pair correlation functions between 

the different spatial components of spins Cff = Cff = g^j ) and Cff = d ^ ^ 

at centers 1 and 2, or the z-component of quadrupolar moment C** = d ^ D ^ '■ 



m z = l -{St + SI) = i{2 exp[-/3(J 2 + 2D)} sinh(2/3#) 

+2exp(-(3D) sinh(/3#) cosh^J^)}, (7) 

= ±(SfS% + S?S%) = "^{2 exp(-pD) cosh(/3tf) smh([3J xy ) 

+2exp[p(J z - 2D) 12} smh(PR)/R}}, (8) 

C\l = (S^SI) = ^{2 exp[-/3( J z + 2D)} cosh(2/3if ) - exp[/3( J z - 2D)} 

-exp[(3(J z - 2D)/2}[cosh((3R) + (J z - 2D) smh((3R)/2R]}, (9) 

C{{ = lm) 2 + (SI) 2 ) = i{2 exp[-(3(J z + 2D)} cosh(2(3H) 
+2exp{-(3D) cosh{(3H) cosh{(3J xy ) + exp[/3(J 2 - 2D)] 
+ exp[(3(J z /2 - D)][cosh((3R) + (J z - 2D) sinh((3R)/2R]}. (10) 

It is worthwhile to mention that the quantities listed above can also be calcu- 
lated from the definition of statistical mean value after performing the trans- 
formation of the relevant quantity into the basis of eigenfunctions, i.e. by 
means of: (...) = J2 a (^a\---\^a)e~ l3Xa (summation runs over all eigenstates). 
This alternative approach is even more general, since it enables to calculate 
the ensemble average for any combination of random spin variables and not 
only for those, which are explicitly involved by the different interaction pa- 
rameters in the total spin Hamiltonian. As an example, we have evaluated 
within this approach the transverse component of quadrupolar moment, the 
quantity, which cannot be calculated by differentiating the Gibbs free energy: 
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dt = ((S*) 2 } = ((Sf) 2 } = -{exp[-/?( J z + 2D)} cosh(2/3H) 
+ 3 exp(-/3L>) cosh(/3#) cosh((3J xy ) + exp[/3(J 2 - 2D)}/2 
+ exp[(3(J z /2 - D)][3cosh(/3i2)/2 - (J z /2 - D)smh((3R)/2R]}. (11) 

It is noticeable that the calculated correlation functions and spatial compo- 
nents of quadrupolar moment provide, in addition to the better understanding 
of the spin order, useful connection to the macroscopically measurable quan- 
tities such as susceptibility, entropy, specific heat, etc. In fact, by exploiting 
the fluctuation-dissipation theorem one finds following relations between the 
spatial components of molar susceptibility and respectively, magnetization, 
correlation functions and quadrupolar moment: 



X:i 



Xz 



■■N. 



Knl 



■■N. 



A f^(Cr i +C--2m 



;i2) 



With regard to earlier results (7)-(ll), the set of Eqs. (12) yields the closed- 
form expressions allowing the exact calculation of the longitudinal (x z ) as 
well as transverse (xx, Xy) components of the molar susceptibility (A^a de- 
notes Avogadro's number). Of course, similar relations can also be established 
between the internal energy, entropy or specific heat on the one side and the 
magnetization, correlation functions and quadrupolar moment on the other 
side. However, these thermodynamic quantities can be readily obtained from 
the Gibbs free energy with the help of standard thermodynamic relations and 
it is therefore meaningless to use such an approach for them. 



3 Results and discussion 



Before proceeding to the discussion of the most interesting analytical and 
numerical results, it should be emphasized that the exact results derived in 
the previous section are rather general, in fact, they hold for any combination 
of intra-dimer coupling constants J xy and J z (both of them can be F or AF), 
arbitrary uniaxial anisotropy D and external magnetic field B. For simplicity, 
we will suppose that both coupling constants are of equal sign, i.e. they are 
either ferromagnetic (J xy , J z < 0), or antiferromagnetic (J xy , J z > 0). 
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3.1 Ground-state behavior 

At first, let us focus on the ground-state properties of spin-1 dimer. Since the 
ground state is definitely determined by the lowest-energy eigenstate, it is just 
a problem of finding the lowest energy level among all possible eigenenergies. 

3.1.1 Zero-field case 

In the zero-field case we deal with a rather simple situation. It can be eas- 
ily proved that the |^o,o) is always the lowest-energy eigenstate when both 
coupling constants are AF (J xy , J z > 0). In the consequence of that, the spin 
arrangement of the spin-1 AF dimer can be characterized with the aid of: 



It can be easily understood that the strength of AF spin alignment within 
different spatial directions depends basically on a competitive influence of J z , 
J xy and D interaction parameters. As could be expected, the easy-axis single- 
ion anisotropy forces spins to align anti-parallel with respect to each other, 
letting D — > — oo indeed gives Cff = — Cff = —1 similarly as it is in the 
case of spin-1 Ising dimer, i.e. under the assumption of the infinitely strong 
easy-axis exchange anisotropy (J z — > oo). On the contrary, there is an obvious 
distinction between the effect of exchange and single-ion anisotropies in the 
easy-plane regime. The easy-plane single-ion anisotropy clearly tries to tend 
both spins in the xy-plane. Apart from Cff = C\{ = retrieved in the limit 
D — > +oo, the |\l/o,o) eigenvector can serve in evidence of this statement in 
that it reduces towards a single |0, 0) state. On the other hand, the population 
of the 'intrinsic AF' | ± 1, =)=1) states does not completely vanish owing to the 
easy-plane exchange anisotropy even if J xy — > oo. Actually, the spin-1 dimers 
are constrained to remain in the superposition of spin states even under this 
condition and the 1 0, 0) state is occupied with the same probability as both 
'intrinsic AF' | ± 1, =)=1) states together. An independent check of this scenario 
also provides C*% = — C\\ = —0.5 retrieved in the J xy — > oo limit. 

Surprisingly, the situation becomes a bit more involved when assuming a fer- 
romagnetic version of the spin-1 dimer (J xy , J z < 0). The doubly degenerate 
energy level that corresponds to both eigenvectors |^2,±2) with a parallel spin 
alignment (spins are oriented either upwards or downwards) is, namely, the 
lowest-lying only in a certain subspace of interaction parameters. Really, the 




(13) 
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easy-plane anisotropy may induce a level crossing between the doublet of F 
eigenstates |*2,±2} and respectively, the eigenstate |*2,o} which can also be 
characterized by set of Eqs. (13). It can be realized that the considered level 
crossing occurs on behalf of the competition between J z , D and J xy when: 

-J l :=\J l z C \=D + ^D* + .P xy1 (for J z < only). (14) 

The lowest-energy level is the doublet of ^,±2} states if | J z \ > \J l z c \ is satisfied, 
otherwise the non-degenerate |*2,o} state possesses the lowest energy. As far 
as the single-ion (exchange) anisotropy is concerned by itself only, the value 
matching the level crossing being -J l z c = 2D (J l z c = J xy ). Accordingly, it can 
be easily understood that an appearance of the 'symmetric' AF eigenstate 
|*2,o} in the region of F coupling constants ( J xy , J z < 0) is strongly reinforced 
when the easy-plane single-ion anisotropy acts simultaneously with the easy- 
plane exchange anisotropy. However, it should be also mentionedthat the |*2,o} 
state is energetically favored before the doublet of F states |^/ 2 ,± 2 ) provided 
that the easy-axis exchange anisotropy is accompanied with much stronger 
easy-plane single-ion anisotropy or naturally, in the reverse case. 

3.1.2 Effect of external field - the magnetization process 

In this part, we will explore the effect of external magnetic field on the 
character of ground-state spin arrangement. To simplify further manipula- 
tions, we will rescale all interaction parameters with respect to the exchange 
coupling \J xy \ by introducing the following set of normalized parameters: 
A = J z /\ J xy \; d = D/\J xy \; h = H/\J xy \ = gzfi^B /\J xy \; t = k^T/\J xy \. 

For illustration, Fig. 1 displays the regions corresponding to the lowest-energy 
eigenstates in the d — h plane for several values of the exchange anisotropy A. 
Upper (lower) panel fits the situation by considering the spin-1 AF (F) dimer; 
the left, central and right panels serve for a comparison of the easy-plane, 
isotropic and easy-axis interaction regimes. Apart from the |*o,o} an d |*2,o} 
eigenstates, which have already been described by means of Eqs. (13), the 
other allowable zero-point eigenstates can be characterized as: 

: = -0.5, Clt = 0.75, C\l = 0.0, C{{ = 0.5, m z = 0.5; (15) 
|tf 2>1 ) : C™ = 0.5, CH = 0.75, C\l = 0.0, C{[ = 0.5, m z = 0.5; (16) 
^2,2) : C% = 0.0, Cft = 0.5, C\l = 1.0, C{{ = 1.0, m z = 1.0. (17) 

Further, the boundaries separating the lowest-energy eigenstates are given by: 
|*o,o>(|*2,o» - |*i,i)(|*2,i» : h = X/2-l + y/(d - A/2)2 + 2; (18) 
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|*i,i)(l*2,i))- |* 2)2 ): h = X + d+l; (19) 
|*o,o)(|*2,o)) - 1*2,2) : h = 3A/4 - d/2 + y/(d - X/2f + 2/2; (20) 

and for completeness, we also write in a parametric form the coordinates of 
special points at which all three boundary lines cross each other: 

A>-2; d t = -(A + l)/(A + 2); h t — (X + 1) 2 /(A + 2). (21) 

Apparently, the zero-temperature magnetization curves are obtained from Fig. 
1 when moving along the vertical fa-axis. If doing so, one gains a stepwise 
magnetization curve with the abrupt change (s) of magnetization at crossing 
field(s) given by Eqs. (18)-(20). It is quite obvious that there are at best two 
distinct types of magnetization curves for the spin-1 AF dimer. Firstly, one 
encounters a single-step magnetization curve with a direct field-induced tran- 
sition from the |*o,o) state to the 1*2,2) state with saturated magnetization. 
Secondly, there also appears a two-step magnetization curve with an interme- 
diate plateau at half of the saturation magnetization in between the above 
mentioned |*o,o) and 1*2,2) states. Evidently, a presence of the intermediate 
plateau reflects the energetical stability of the |*i,i) state in the range of 
moderate fields. As could be expected, the easy-axis (easy-plane) anisotropy 
prefers the single-step (two-step) magnetization curves no matter whether it 
is caused by the single-ion anisotropy or exchange anisotropy term. Finally, it 
is noteworthy that the purely isotropic spin-1 AF dimer (A = 1, d — 0) shows 
the double-plateau magnetization curve. 

As far as the spin-1 F dimer is concerned, the stepwise magnetization curves 
appear just as the claim | J z \ < \J l z c \ is fulfilled because of the strong enough 
easy-plane anisotropy. Despite the magnetization curves seem in such a case 
nearly the same as those of the spin-1 AF dimer, there are nevertheless some 
principal differences. At first, the zero-field as well as intermediate-plateau 
states do not occur owing to a presence of the |*o,o) and |*i,i) states, but 
on behalf of the |*2,o) and |*2,i) states. Though these eigenfunctions differ 
from each other just in the symmetry, this change gives rise to a change of the 
character of spin correlations in the xy-pleme (see Eqs. (13), (15) and (16)). 
Another important difference emerges when considering the combination of 
the easy-axis exchange anisotropy and easy-plane single-ion anisotropy. As 
the easy-axis exchange anisotropy strengthens, the boundary lines separating 
the intermediate-plateau state |*2,i) merge together until the plateau state 
wholly disappears from the magnetization curve below A < —2. This behavior 
is in a sharp contrast with what is observed in the spin-1 AF dimer, where the 
plateau state |*i,i) does not vanish due to the easy-plane single- ion anisotropy 
nothwithstanding of the effect of easy-axis exchange anisotropy. 

To conclude our analysis of the ground state, we depict in Fig. 2 the Zeeman 
splitting of energy levels in the spin-1 AF dimer by selecting several combina- 
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tions of the anisotropy parameters d and A. Fig. 2a shows a well-known picture 
matching the situation in the isotropic spin-1 dimer; the energy spectrum con- 
sists of the singlet (|^ ,o>), triplet (|*i )0 ), |*i,±i» and quintet (|* 2) o), |*2,±i), 
I ^2, ±2)) energy levels and the level degeneracy is entirely lifted by any non- 
zero external field. To provide a deeper insight, Fig. 2b (Fig. 2c) illustrates 
the effect of easy-plane (easy-axis) single-ion anisotropy and Fig. 2d (Fig. 2e) 
the effect of easy-plane (easy-axis) exchange anisotropy on the overall trends 
in the energy spectrum. As one can see, both anisotropy parameters are re- 
sponsible for the zero-field splitting: the triplet level is consequently splitted 
into the doublet |^i,±i) and the singlet |^i,o) level, whereas the quintet level 
splits into two doublets ^,±1), ^,±2) and one singlet |\&2,o) level. Notice 
that the zero-field splitting induced by the exchange anisotropy and the one 
originating from the single-ion anisotropy differ from each other even if both 
anisotropies are assumed to be of the same easy-axis (easy-plane) type. Appar- 
ently, the aforementioned difference rests in the splitting of the quintet level, 
while the triplet level is, on the contrary, splitted in the qualitatively same 
manner. Finally, it should be pointed out that by reverting the sign of vertical 
(energy) axis, one obtains the energy spectrum for the spin-1 F dimer with 
the interaction parameters given by the transformation (X,d) — > (—A, —d). 



3.2 Finite-temperature behavior 

Now, we shall pay our attention to the finite-temperature behavior of the 
magnetization and susceptibility, the quantities, which characterize a response 
of the magnetic system with respect to the external field. Both possible kinds of 
magnetization curves (single-step as well as two-step) are illustrated in Fig. 3. 
It is necessary to emphasize that there are no real jumps in the magnetization 
vs. field dependence at any finite temperature and that the true discontinuities 
occur merely at zero temperature as the external field reaches some of its level- 
crossing values (18)-(20). At finite but sufficiently low temperatures, however, 
an extremely steep but continuous increase of the magnetization observable 
in the vicinity of the level-crossing fields strongly resembles the ground-state 
magnetization jumps. It is easy to understand from Fig. 3 that the abrupt 
change of magnetization is gradually smeared out by the thermal fluctuations 
as the temperature raises and in the consequence of that, the magnetization 
curve acquires at sufficiently high temperatures a typical magnetization profile 
of the Langevin-Brillouin type. It is worthy to mention that the depicted 
magnetization curves render the temperature effect on the overall trends in 
the magnetization profile even if A 7^ 1. 

Furthermore, the temperature variations of the longitudinal (x z ) and trans- 
verse (xx) susceptibilities are displayed in Fig. 4 for several combinations of 
the anisotropy parameters A and d. As one can see from Fig. 4a, a round 
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maximum in the t — Xz dependence flattens and shifts towards higher tem- 
peratures when the easy-plane single-ion anisotropy strengthens. Conversely, 
the height of maximum is almost independent of the easy-axis single-ion 
anisotropy strength and there is only a slight shift of the maximum posi- 
tion by changing this type of magnetic anisotropy (Fig. 4b). As far as the 
thermal variation of Xx is concerned, the transverse susceptibility diverges by 
approaching the zero temperature whenever d 7^ and/or A ^ 1, i.e. whenever 
there is at least small magnetic anisotropy. It is quite apparent from Figs. 4c- 
d that there is a quantitative rather than qualitative difference between the 
effect of easy-plane and easy-axis single-ion anisotropies on the thermal de- 
pendence of transverse susceptibility. At relatively weak single-ion anisotropies 
(\d\ < 0.75), the t — Xx dependence exhibits, in addition to the familiar round 
maximum, the successive minimum before showing the final low-temperature 
divergence under the temperature suppression. The round minimum in the 
t — Xx dependence gradually flattens as the single-ion anisotropy strengthens 
and the transverse susceptibility becomes a monotonical decreasing function 
of temperature above \d\ > 0.75. Finally, we just quote that the addition of 
the exchange anisotropy has a very similar effect on the thermal dependences 
of longitudinal and transverse susceptibilities as the single-ion anisotropy. 



4 Magnetic behavior of homodinuclear nickel(II) complex 
[Ni 2 (Medpt) 2 (^-ox)(H 2 0) 2 ](C10 4 ) 2 .2H 2 (NAOC) 

In this part, the obtained theoretical results will be confronted with the rel- 
evant experimental measurements performed on the single-crystal sample of 
[Ni 2 (Medpt) 2 (/i-ox)(H 2 0) 2 ](C10 4 ) 2 .2H 2 (hereafter abbreviated as NAOC), 
which is regarded as a typical representative of the spin-1 dimeric compound 
[19]. A schematic view on the dinuclear unit, relative positions of non-ligating 
perchlorate counter-anion and water molecules are presented in Fig. 5. Both 
octahedrally coordinated nickel(II) atoms are linked via bis-chelating oxalato 
group and the rest of their coordination sphere is occupied by the block- 
ing tridentate amine (Medpt) and ligating water molecule. Accordingly, the 
single-crystal sample of NAOC consists of an assembly of discrete dinuclear 
entities, which are connected through the hydrogen bonds and/or van der 
Waals forces to the molecular crystal. It can be viewed from Fig. 5 that the 
NAOC molecule is centrosymmetric with an inversion center at the midpoint 
of the bridging oxalato group. Thereupon, the antisymmetric Dzyaloshinski- 
Moriya interaction [20] (a first-order effect of spin-orbit coupling and exchange 
interaction) should be entirely neglected as a possible source of the magnetic 
anisotropy in the NAOC. The exchange anisotropy, i.e. the anisotropy in a 
symmetric (pseudodipolar) exchange interaction, could be then anticipated to 
be the third strongest source of magnetic anisotropy beyond the single-ion 
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anisotropy and the g-factor anisotropy. 

Before proceeding to a comparison of the theoretical and experimental results 
it should be noted here that the magnetic measurements performed on the 
NAOC sample have been primarily reported by some of the present authors 
in the earlier publications [21]- [22] to which the interested reader is referred 
for a closer experimental details. Notice that both magnetization as well as 
susceptibility data reported on previously imply strongly anisotropic magnetic 
behavior [21] and the marked anisotropy manifests itself also in the peculiar 
angular dependence of level-crossing fields (see Fig. 1 of Ref. [22]). By taking 
into account the results of Ref. [22], we choose the crystallographic c*-axis 
as the principal axis because c* being the easy magnetization axis. The crys- 
tallographic a-, b- and c*-axis are consequently defined to be the x-, y- and 
z-axis of the effective spin Hamiltonian. Finally, we should also remark that 
the high-field magnetization data used in our further analysis are identical 
with those reported on earlier [21], while the original susceptibility data were 
found erroneous and we report here a revised set of data for them. 

Magnetization normalized with respect to its saturation value are plotted in 
Fig. 6 against the external magnetic field applied along the c*- and a-axis. 
Unfortunately, there does not exist an exact analytical expression for the mag- 
netization of spin-1 dimer when the external field is applied perpendicular to 
the principal axis. Adopting the least-square fitting based upon a minimali- 
sation of the function R = Y.i.^T 1 ~ A° bs ) 2 / *E(A° bs ) 2 , both displayed magne- 
tization curves were therefore fitted through the relation (6) derived for the 
magnetization along the principal direction. At first, we decided to neglect the 
exchange anisotropy by keeping both exchange parameters equal to each other 
(Jxy — Jz — J) and only the gyro magnetic ratio g together with the single-ion 
anisotropy parameter D were allowed to vary within different spatial direc- 
tions. Under these conditions, the best simultaneous fit of all magnetization 
data is obtained for the following set of parameters (Fig. 6): J/k B = 30.66K 
and D c */k B = -12.48K, g c * = 2.28 (R c * = 0.0019) for c*-axis and respectively, 
Da/ks = 4.91K, g a = 2.24 (R a = 0.0067) for a-axis. In agreement with our 
expectations, the magnetization fit is substantially better for the c*-axis than 
for the a-axis (R c * < R a ) and the anisotropy parameter D has an opposite 
sign for the easy (c*) and hard (a) magnetization axis. 

On account of a relatively strong easy-axis single-ion anisotropy, the interme- 
diate plateau shrinks to a very narrow field range when the external field is 
applied along the easy magnetization axis (Fig. 6a). Moreover, it is interesting 
to note that the obtained single-ion anisotropy strength D c * / J —0.41 is 
rather close to the boundary value D t /J = — | below which the intermediate 
plateau should completely disappear from the magnetization curve according 
to the theoretical prediction based on Eq. (21). The inverse effect is observed, 
however, when applying the external field normal to the easy magnetization 
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axis (Fig. 6b). In this case, the field range corresponding to the intermediate 
plateau largely extends not only owing to a much higher saturation field, but 
also due to a slightly lower field required for attaining the plateau state (see 
the insert of Fig. 6b). Finally, it is worthy to note that we found a strong evi- 
dence that the magnetic behavior observed in the NAOC can be undoubtedly 
attributed purely to the single- ion anisotropy effect. The exchange anisotropy 
is unlikely to significantly contamine the magnetochemistry of the NAOC, 
actually, the improvement of R a and R* quality factors accomplished by the 
inclusion of this anisotropic term does not exceed 0.5%. 

The temperature dependences of susceptibility measured along the a-, b- and 
c*-axis are compared in Fig. 7 with the relevant theoretical prediction for 
the longitudinal and transverse susceptibility. Notice that the experimental 
susceptibilities were corrected for the diamagnetic contribution according to 
Pascal tables. Under the restriction to the isotropic intra-dimer interaction J, 
the best simultaneous fit of all susceptibility data was obtained for this unique 
set of interaction parameters: J/k-Q = 30.89K, D/k-Q = — 5.61K, g a = 2.34 
(R a = 0.0028), g b = 2.30 (R b = 0.0019) and g c * = 2.37 (R c * = 0.0050). 
Evidently, there is a rather reasonable accord between the values of J, g a , g b 
and g c * parameters listed in the fitting set of magnetization and susceptibility, 
as a matter of fact, the latter values are only slightly enhanced with respect 
to the former ones. The most significant discrepancy thus rest in the strength 
of uniaxial single-ion anisotropy D, which is more than twice weaker for the 
susceptibility fit than for the magnetization fit. 

There are strong structural indications that the inconsistency in determining 
the uniaxial zero-field splitting parameter D can be related to higher-order 
anisotropy terms, which were entirely neglected in our study. The ligand asym- 
metry in the equatorial plane of Ni 2+ ions (see Fig. 5) indeed implies a possi- 
bly non-negligible biaxial single-ion anisotropy and the numerical calculations 
involving this anisotropic term strongly support this concept [22]. Unfortu- 
nately, the influence of biaxial anisotropy was out of scope of the present study, 
since this term prevents an exact analytical treatment used here. Further, it 
turns out that the monotonical decrease of susceptibility observed upon cool- 
ing ceases below 5K independently of the spatial direction, what indicates a 
small fraction of non-coupled nickel (II) impurities to be present in the sample. 
However, the best simultaneous fit of all susceptibility data which accounts 
also for the paramagnetic impurities considerably reduces just an error in 
the quality factor R c * (from 0.0050 down to 0.0019), while the other quality 
factors and fitting parameters remain by the estimated concentration of impu- 
rities [p ~ 0.34%) almost unaltered. The most substantial deviations between 
the theoretical and experimental curves thus appear in the high-temperature 
region, where the theoretical curves lie slightly below the experimental ones. 
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5 Conclusion 



In the present article, the particular attention has been devoted to a compara- 
tive study of the uniaxial single- ion anisotropy and exchange anisotropy as pos- 
sible sources of the magnetic anisotropy in the spin-1 dimer. Exact analytical 
results for the Gibbs free energy, longitudinal magnetization, correlation func- 
tions, quadrupolar moment, longitudinal and transverse susceptibility were 
derived and on the basis of these results, a simple derivation of other quantities 
such as internal energy, entropy, or specific heat, can easily be accomplished. 
It is worthy to notice that all aforementioned results can be rather straightfor- 
wardly extended to account for an effect of the inter-molecular (inter-dimer) 
interaction treated within the mean-field approximation. 

The most important result stemming from our study consists in the derivation 
of exact closed-form expressions for the longitudinal and transverse suscepti- 
bilities by means of the fluctuation-dissipation theorem, which relates these 
quantities to the magnetization, correlation function and quadrupolar mo- 
ment. To the best of our knowledge, we are not aware of any work providing 
an exact closed-form formula for the transverse susceptibility of the anisotropic 
spin-1 dimer. In addition, we found a convincing evidence that the transverse 
susceptibility of spin-1 dimer diverges as T — * whenever there is an arbitrary 
but non-zero magnetic anisotropy. Finally, the accurate condition allocating 
the 'symmetric' AF ground state of spin-1 F dimer was determined (14) and 
its physical background was discussed in detail. 

As far as the single-ion anisotropy and exchange anisotropy effects are com- 
pared, the field dependence of magnetization as well as the temperature de- 
pendence of susceptibility are influenced by these anisotropic terms in a qual- 
itatively same manner. It is therefore rather intricate to find out, sometimes, 
whether the overall magnetic anisotropy is attributable purely to one source, 
or arises as a sum of more contributions. It should be mentioned here, how- 
ever, the high-field ESR measurement provides a very efficient alternative tool 
that enables to distinguish between the different magnetic anisotropies, since 
it reflects the overall energy spectrum depending qualitatively on the type of 
magnetic anisotropy as one can understand from the set of eigenenergies (4). 

The obtained theoretical results were also compared with the available exper- 
imental data of the NAOC, which is regarded to be a typical representative 
of the spin-1 AF dimer compound. There is a strong evidence that the zero- 
field splitting induced by the uniaxial single-ion anisotropy is by far the most 
dominating phenomenon which determines the overall magnetic anisotropy 
in the NAOC compound. Nevertheless, it is interesting to mention that the 
strength of uniaxial single-ion anisotropy evaluated from the magnetization 
fit is quite close to the boundary value below which the system should ex- 
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hibit a striking spatially-dependent magnetization process. Subsequently, one 
should encounter the single-step magnetization curve when applying the exter- 
nal field along the easy magnetization axis, while the two-step magnetization 
curve should be expected to occur when the field is applied perpendiculary to 
it. Although the design of molecule-based magnetic materials with a tunable 
strength of the single- ion anisotropy is far from a routine target at present, a 
large amount of oxalato-bridged dinuclear nickel complexes [6] supports our 
hope that this family of compounds will provide a suitable candidate for dis- 
playing such an interesting phenomenon. 
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Figure captions 

Fig. 1 The lowest-energy eigenstates in the d — h plane for several values of ex- 
change anisotropy A. Upper (lower) panels fits the situation in the spin-1 
AF (F) dimer, the left, central and right panels serve for a comparison of 
the easy-plane, isotropic and easy-axis interaction regimes, respectively. 

Fig. 2 Zeeman's splitting of the complete energy spectrum shown for several com- 
binations of the anisotropy parameters A and d: a) A = 1.0, d = 0.0; b) 
A = 1.0, d = 0.5; c) A = 1.0, d = -0.5; d) A = 0.5, d = 0.0; and e) 
A = 2.0, d = 0.0. For clarity, the indices of eigenvectors are shown only. 

Fig. 3 The typical examples of two-step and single-step magnetization curves dis- 
played at various dimensionless temperatures t = 0.00,0.05,0.10,0.20,0.50 
in ascending order along the direction of arrows. The parameters used are 
A = 1.0, d = 0.50 (Fig. 3a) and A = 1.0, d = -0.75 (Fig. 3b), the magneti- 
zation is scaled in the g z [iB units. 

Fig. 4 The thermal dependence of the zero-field susceptibility scaled in the N ^(g n B ) 2 / k B 
units when the parameter A is fixed and the parameter d varies. Upper 
(lower) panel shows the longitudinal (transverse) susceptibility, the left 
(right) panels illustrates the situation under the easy-plane (easy-axis) single- 
ion anisotropy. 

Fig. 5 Schematic representation of the NAOC dinuclear unit, non-ligating per- 
chlorate counter anion and water molecules. The NAOC molecule has an 
inversion center at the midpoint of the oxalato bridging group. 

Fig. 6 The high-field magnetization curves measured in the pulsed magnetic field 
along the a- and c*-axis at very low temperature T = 1.3K. The best si- 
multaneous fit, determined as described in the text, was obtained for the 
following set of fitting parameters: J/k B = 30.66-fT and D c */k B = — 12.48-fT, 
g c * = 2.28 (R c * = 0.0019) for c*-axis and respectively, D a /k B = A.91K, 
g a = 2.24 (R a = 0.0067) for a-axis. The solid (dashed) lines show the ex- 
perimental (theoretical) curves, the insert depicts both experimental mag- 
netization curves in the same scale. 

Fig. 7 The thermal variations of magnetic susceptibility measured along the a-, 
b- and c*-axis plotted together with the best simultaneous fit obtained for: 
J/k B = 30.89fT, D/k B = -5MK, g a = 2.34 (R a = 0.0019), g b = 2.30 
(R b = 0.0021) and g c * = 2.37 (R c * = 0.0050). For clarity, the theoretical 
predictions for the longitudinal and transverse susceptibility representing, 
respectively, the susceptibility measured along the c*- and 6-axis are shown 
only. 
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Fig.2 Strecka et al 
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Fig. 5 Strecka et al. 




Fig. 6 Strecka et al. 



